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We present S factor data obtained from the Coulomb dissociation of 83 MeV/nucleon *B, and ana- 
lyze ^Be longitudinal momentum distributions measured at 44 and 81 MeV/nucleon using a potential 
model, first-order perturbation theory, and dynamical solution of the time-dependent Schrodinger 
equation. Comparing our results with independent continuum-discretized coupled channels calcu- 
lations, we study the reaction model and beam energy dependence of the E2 contribution to the 
dissociation cross section. By fitting radiative capture and Coulomb breakup data taken below rela- 
tive energies of 400 keV with potential models constrained by ^Li -I- n and ^Be -I- p elastic scattering 
data, we examine the mutual consistency of recent S17 measurements and obtain a recommended 
value for 817(0) of 18.6 ±0.4 (experimental) ± 1.1 (extrapolation) eV b (Icr). This result is in good 
agreement with recent experimental determinations of the asymptotic normalization coefficient of 
the valence proton wave function in *B. 
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I. INTRODUCTION 

With the discovery that a large fraction of the high- 
energy electron neutrinos emitted in the (i'^ decay of 
in the sun are transformed into other active neutrino fla- 
vors on their way to terrestrial detectors and, consis- 
tent with CPT symmetry, that reactor-produced antineu- 
trinos also oscillate ^ , a robust solution of solar neutrino 
problem seems to be in hand. Although the SNO neutral 
current measurement of the ^B solar neutrino flux is con- 
sistent with the predictions of the standard solar model 
[sl, its precision is stiU far from that required to ascer- 
tain if there is transformation of solar electron neutrinos 
into sterile neutrinos |^ |^ . In order to reach this goal, 
the precision of both the experimental measurement and 
the theoretical prediction must be improved. Through 
the addition of NaCl and "^He proportional counters to 
the SNO detector, its neutral current measurement will 
be made much more precise. In order to improve the 
precision of the theoretical prediction of the *B neutrino 
flux, the rate of the reaction that produces ®B in the sun, 
^Be(p, 7)®B, must be better determined. Currently, the 
uncertainty on the zero-energy astrophysical S factor for 
this reaction, 817(0) makes the largest contribution to 
the theoretical error budget. 

Despite the profusion of recent determinations of 
Si7(0) |SllSSES'Tl','T5,lil|ll|ll|ll|llia, the cur- 
rent recommendation 19] is based on the results of a sin- 
gle experimental measurement |23| ■ Recently, a weighted 
mean of such determinations based on both direct and in- 



direct methods was given [2l|. Here, we present data ob- 
tained via the indirect method of Coulomb dissociation 
p^ . offer theoretical analyses of these and other data, 
and assess the current state of knowledge of 817(0). 

At solar energies, the radiative capture ''Be(p, 7)^B is 
entirely El dominated. The Coulomb dissociation of ^B 
at intermediate beam energies also proceeds predomi- 
nantly via E\ transitions, but the much larger flux of 
E2 virtual photons leads to a finite E2 contribution to 
the dissociation cross section. This complicates efforts 
to extract the El strength that determines the radiative 
capture rate. 8ince the first measurement of the Coulomb 
dissociation of *B j23], the role of E2 transitions in this 
process has been rather controversial. Experimental ef- 
forts to determine the E2 contribution to the breakup 
cross section have focused on the relative energy and ^B 
scattering angle distributions 0, 01 ; the angular distri- 
butions 18] and longitudinal momentum distributions of 
the emitted protons ^^li ^"^^ the longitudinal momen- 
tum distributions of the ^Be fragments |2^. The 
relative energy and ^B scattering angle distributions are 
relatively insensitive to E2 strength, since the El and 
E2 contributions sum incoherently ,27,]. In contrast, the 
fragment longitudinal momentum and angular distribu- 
tions are better suited to this task because interference 
between El and E2 transition amplitudes produces strik- 
ing asymmetries in these distributions |28| . 



II. LONGITUDINAL MOMENTUM 
DISTRIBUTIONS AND E2 STRENGTH 
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A first-order perturbation theory analysis of ^Be longi- 
tudinal momentum distributions found that the E2 ma- 
trix elements implicit in the potential model used had to 
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be scaled by 0.7 in order to best fit the measured distri- 
butions [2l|. In that work, the notion that higher-order 
dynamical effects quench the E2 strength predicted by 
first-order theories was advocated. Support for this 
idea was provided by concurrent |i21j| and subsequent 
continuum-discretized coupled channels (CDCC) calcu- 
lations. The perturbative calculations employ a sim- 
ple, single-particle model for the structure of *B, and 
calculate the electromagnetic dissociation cross section 
assuming single photon exchange, neglecting nuclear- 
induced breakup. The CDCC calculations use a similar 
structure model but include electromagnetic and nuclear 
matrix elements to all orders. In Ref. j^^, the authors 
found that a scaling of the E2 matrix elements by a fac- 
tor of 1.6 was required for the best description of the data 
using CDCC. Hence the E2 strength had to be reduced 
in first-order calculations, and enhanced in all-order cal- 
culations, relative to the original structure model predic- 
tions. 

Here we study the same longitudinal momentum distri- 
butions ^2?\ using two approaches, first-order perturba- 
tion theory and dynamical solution of the time-dependent 
Schrodinger equation ,30] ■ Both methods employ only 
Coulomb matrix elements, and assume the same single- 
particle potential model. In this model, the p-wave po- 
tential depth has been fixed by the ®B binding energy, 
and the depths of the potential for all other partial waves 
have been fixed by the s-wave scattering lengths for chan- 
nel spin 1 and 2 in the isospin mirror ^Li -I- n system (3lj . 
The potential radius is 2.5 fm, the diffuseness 0.65 fm, 
and the p-wave potential depth is 43.183 MeV. For s, d, 
and /-waves, the potential depths are 43.857 and 52.597 
MeV for S — 1 and 5 = 2 respectively. Spectroscopic 
factors for the two spin configurations were taken from 
Ref. after a centcr-of-mass correction, and are 0.3231 
for 5 = 1 and 0.8572 for S — 2. These theoretical spec- 
troscopic factors agree very well with other shell model 
calculations and cluster model calculations (33| . Our 
potential model and that used in Ref. 29] differ from 
that of Ref. 2^ in that they lack a spin-orbit interac- 
tion. Accordingly, our potential model includes no 1+ 
resonance at 0.6 MeV and lacks resonant E2 strength in 
this region. The present model gives S'_e2/5'_ei(0.6 MeV) 
= 5.9 X 10^"*, compared with the 9.5 x 10"'* predicted by 
the model of Ref. 28] . These two values span the range 
of theoretical predictions made using disparate methods 
cataloged in Ref. [2l| . 

We have scaled the E2 matrix elements calculated with 
our potential model, and then performed perturbative 
and dynamical calculations with the different E2 scaling 
factors. In the dynamical calculations the wave function 
of the ''Be -I- p system was expanded into partial waves, 
taking relative orbital angular momenta £ < S into ac- 
count. The radial wave function for each partial wave, 
specified by £ and its projection m, was discretized on 
a mesh with exponentially increasing step size covering 
radii between and 900 fm. Starting with the ground 
state of ^B, the evolution of the coupled radial wave 



functions was followed by applying a unitary approxima- 
tion of the time-evolution operator in steps of 5 fm/c. A 
smooth switch-on and switch-off of the Coulomb poten- 
tial for large distances was introduced to avoid spurious 
excitations. Only multipoles A = 1 and 2 were considered 
in the expansion of the time-dependent Coulomb inter- 
action while the *B center of mass followed a hyperbolic 
Coulomb trajectory with respect to the target. The final 
time-evolved wave function was projected onto the rele- 
vant continuum wave functions to obtain the excitation 
amplitudes and the cross sections. Allowing the overall 
normalization to vary freely and convoluting the pertur- 
bative and dynamical calculations with the experimental 
momentum resolution of 5 MeV/c, we have determined 
the E2 scaling factor that minimizes the value for the 
central six 44 MeV/nucleon data points and the central 
five 81 MeV/nucleon data points of the measured longi- 
tudinal momentum distributions. 

Fig- n shows the results of the best- fit calculations for 
the breakup of ^B on Pb at 44 MeV/nucleon, while Fig. 
121 shows the 81 MeV/nucleon results. Both the perturba- 
tive and the dynamical calculations satisfactorily repro- 
duce the measured asymmetries in the central regions of 
the longitudinal momentum distributions, though neither 
calculation describes the widths of the 81 MeV/nucleon 
distributions very well; the measured distributions are 
wider. As the tails of the momentum distributions cor- 
respond to high *B excitation energies, it is possible that 
our neglect of the nuclear interaction, which is relatively 
more important at high excitation energies than low, 
is responsible for this discrepancy. The required over- 
all normalization factors range from 0.7 to 0.95 for the 
various angle cuts, beam energies, and reaction models. 
Consistent with earlier findings [2ll 123, 123, ya] , higher- 
order effects present in the dynamical calculations but ne- 
glected in the first-order perturbative calculations tend to 
reduce the asymmetry predicted for a given E2 strength. 
Hence smaller E2 scaling factors are required for the 
perturbative calculations than for the dynamical calcu- 
lations. In addition, there is a hint of a beam energy de- 
pendence apart from that inherent in standard Coulomb 
excitation theory, with less E2 strength required to ex- 
plain the asymmetries observed at the higher beam en- 
ergy than at the lower. Such an effect, which has thus 
far eluded theoretical explanation, may be responsible 
for the reported lack of evidence for E2 contributions to 
Coulomb breakup at 254 MeV/nucleon [8, 18]. 

The best-fit E2 matrix element scaling factors and 
their la uncertainties obtained using the two reaction 
models are shown in Tabled It is of interest to compare 
the results found here with those reported in Ref. |2]| . 
which were deduced from a perturbative analysis of the 
3.5° 44 MeV/nucleon longitudinal momentum distribu- 
tion. The best-fit scaling factor obtained from compar- 
ison of our perturbative calculation with the same data 
here implies Se2/Sei{0-6 MeV) = 5.9^2:1 x lO""^, which 
is consistent with the value of 4.71^ 3 ^ 10""' given in Ref. 
[21I I. From this we conclude that the results of pertur- 
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FIG. 1: Measured '^Be longitudinal momentum distributions 
from the Coulomb dissociation of 44 MeV/nucleon on Pb 
with '^Be scattering angle cuts of 1.5°, 2.4°, and 3.5°. Only 
relative errors are shown. The solid curves in panel (a) are the 
results of first-order perturbation theory calculations, while 
panel (b) shows dynamical calculations. In both cases the E2 
matrix elements have been scaled by the indicated factor. 
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FIG. 2: Same as Fig. [T] but for 81 MeV/nucleon *B on Pb 
with ^Be scattering angle cuts of 1.0°, 1.5°, 2.0°, and 2.5°. 



bative analyses favored thus far by workers in the field 
are relatively insensitive to the details of the underlying 
structure model. This analysis corroborates the case for 
including El contributions in theoretical descriptions of 
the Coulomb breakup of at beam energies below 100 
MeV/nucleon. 

One can also compare the results obtained here with 
dynamical calculations to those found using CDCC cal- 
culations in Ref. [2^ . Those authors studied both the 44 
MeV/nucleon and 81 MeV/nucleon longitudinal momen- 
tum distributions. Although no uncertainty was spec- 
ified, the best-fit scaling value reported in Ref. |23| is 



TABLE I: Best-fit scaling factors of E2 matrix elements pre- 
dicted by our potential model obtained from fitting the lon- 
gitudinal momentum distributions of '^Be fragments from the 
Coulomb dissociation of *B on Pb. The uncertainties are la 
values. 



Beam Energy (MeV/nucleon) 


Calculation 


Scaling Factor 


44 


perturbative 


1.01 ±0.21 


44 


dynamical 


1.6 ± 0.2 


81 


perturbative 


0.63 ±0.12 


81 


dynamical 


1.2 ± 0.2 



1.6, in perfect agreement with the value of 1.6 ± 0.2 
found from the dynamical analysis of 44 MeV/nucleon 
data here. Our result for the 81 MeV/nucleon data is 
somewhat smaller, but if the CDCC scaling factor uncer- 
tainty is comparable to that found here, the two results 
agree reasonably well at this beam energy also. More- 
over, the fact that nuclear matrix elements are included 
in the CDCC calculations but not in our dynamical cal- 
culations suggests that nuclear-induced breakup neither 
significantly enhances nor reduces the asymmetries of the 
Coulomb dissociation calculations for the large impact 
parameters considered here. Since the potential model 
assumed in the CDCC calculations is very similar to that 
employed here, we conclude that these two very differ- 
ent methods predict the same infiuence of higher-order 
effects. We therefore consider these predictions to be ro- 
bust, and conclude that the practice of scaling the E2 
matrix elements in first-order perturbative calculations 
is a reasonable and practical way to proceed with the 
analysis of '^B Coulomb dissociation experiments. 



III. DISSOCIATION CROSS SECTION AND S17 

In Ref. ^2 1 9.n exclusive measurement of the Coulomb 
dissociation of 83 MeV/u ^B on a Pb target was described 
and experimental data were presented in graphical form, 
including the Coulomb dissociation cross section for ^B 
scattering angles < 1.77°. This experimentally measured 
cross section was interpreted in the context of first-order 
perturbation theory in order to obtain 817(0). Here we 
give the energy-dependent S factor derived from these 
data in both graphical and tabular form. Table ^ con- 
tains the measured cross sections and their total Icr un- 
certainties. The systematic uncertainty common to each 
point amounts to 7.1%, and includes contributions from 
the target thickness (1%), beam intensity (2.6%), mo- 
mentum calibration (4.2%), and the size of the E2 com- 
ponent (5%). 

Astrophysical S factors derived from the experimen- 
tal cross section data are shown in Fig. 0] and Table ^1 
These S factors have been corrected for the E2 contri- 
bution to the breakup cross section in the following way. 
The measured cross section was multiplied by the frac- 



TABLE 11: Experimentally measured cross section for the 
Coulomb dissociation of 83 MeV/u *B on Pb for *B labo- 
ratory scattering angles of 1.77° and less (impact parameters 
> 30 fm), and astrophysical S factor for the '^Be(p, 7)*B re- 
action. Uncertainties are la values, and the systematic un- 
certainty common to each point is 7.1%. The relative energy 
bins over which the cross sections and S factors have been 
averaged extend to the midpoints of the intervals separating 
adjacent data points. 



Erei (MeV) da/dErei (mb/McV) Si7(E^ei) (eV b) 

0.064 9.0 ± 2.8 

0.192 50.4 ± 4.9 15.2 ± 1.5 

0.384 113.2 ± 10.1 16.1 ± 1.4 

0.768 116.1 ± 10.9 19.3 ± 1.8 

1.280 61.2 ± 6.8 23.3 ± 2.6 

1.792 41.9 ± 5.9 35.2 ± 5.0 



tion of the total electromagnetic dissociation cross section 
attributable to El transitions to obtain the El cross sec- 
tion. This El fraction was calculated using the Esbensen 
and Bertsch structure model ^231 and first-order pertur- 
bation theory, calibrated by comparison with the mea- 
sured longitudinal momentum distributions as described 
in Ref. [21|. Ml transitions were included by folding 
the S factor measured at the 1"*" resonance by Filippone 
et al. [231 with the calculated virtual photon spectrum. 
The energy- dependent S factor was deduced from the El 
cross section u sing semiclassical first-order Coulomb exci- 
tation theory '3d|, taking into account the experimental 
relative energy resolution. A comparison of theoretical 
predictions of E2 strength found that the smallest 
prediction differed from the largest by less than 40%. Re- 
cent work that examined the influence of structure model 
assumptions on the calculated El fraction [sJl is consis- 
tent with this finding. Since E2 transitions account for 
only 5% of the measured cross section in the perturba- 
tive analysis applied to the data of Ref. the 5% un- 
certainty attributed here to the E2 component is quite 
conservative. Owing to our own findings regarding the 
size of higher-order dynamical effects at very low rela- 
tive energies as well as those of others [s^, we have not 
extracted an S factor from the lowest energy data point. 

In order to minimize the uncertainties due to nuclear 
structure in extrapolating the astrophysical S factor to 
zero energy, it has been recommended that only data ob- 
tained below 400 keV be considered 's^ . This conclusion 
is supported by our own potential model calculations, 
which are shown in Fig. In addition to the potential 
model calculation described above, we have performed 
calculations with the same spectroscopic factors but with 
the s, d, and /-wave potential depths adjusted to re- 
produce the scattering lengths measured in the ^Be -I- p 
system ji^. The potential depths for the central val- 
ues of the measured scattering lengths in this case are 
24.283 MeV and 51.898 MeV for = 1 and 5 = 2 respec- 
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FIG. 3: (Color online) Calculated astrophysical S factors 
for the ^Be(p, 7)*B reaction based on potential models con- 
strained by ^Be -I- p and ^Li -I- n elastic scattering data. The 
1(7 upper and lower limits on the measured scattering lengths 
in the ^Be -I- p system were used to construct parameter sets 
that yield the curves labeled as upper and lower limits. The 
corresponding curves for the ^Li -|- n case are indistinguish- 
able from that based on the central values of these more pre- 
cisely measured scattering lengths. 



tively. The scattering lengths for the ^Be -I- p system are 
not as precisely known as those for the mirror ''Li -I- n 
system; the measured ^Be -I- p and ''Li -|- n scattering 
lengths for the dominant 5 = 2 channel are consistent 
with isospin symmetry, but there is an inconsistency in 
the S — 1 channel that requires further study. Since 
the ''Li -|- n scattering lengths have been measured quite 
precisely, their uncertainties do not significantly affect 
the shape of the S factor calculated with these poten- 
tial model parameters. However, the ^Be -I- p scattering 
lengths have large experimental uncertainties which im- 
ply a correspondingly large uncertainty in the predicted 
shape of the S factor when using a potential model tuned 
to reproduce these elastic scattering data. In Fig. O we 
show the results of the potential model calculations with 
parameter choices that match the measured scattering 
lengths in the two systems. For the ^Be -I- p system, 
we show the potential calculations corresponding to the 
central value and Icr upper and lower limits on the mea- 
sured scattering lengths, while for the ^Li + n system we 
only show the central value calculation because the up- 
per and lower limit curves are indistinguishable from it. 
The extrapolation uncertainty due to the calculated spec- 
troscopic factors is insignificant compared with that due 
to the uncertainty in potential model parameter choices 
corresponding to the different scattering lengths. 

In Fig. ^ the two potential model calculations are 
shown along with the 83 MeV/nucleon Coulomb breakup 
S factor data. The solid line is the potential model with 
^Li -I- n parameters, while the dotted curve is the (cen- 
tral) ^Be -|- p parameter set potential model calculation. 
Both calculations have been scaled to optimally fit the 
data points below 400 keV, and both reproduce the data 
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FIG. 4: (Color online) Astrophysical S factor for the 
^Be(p, 7)*B reaction obtained from the Coulomb dissociation 
of 83 MeV/u *B on Pb. Potential model calculations using 
parameters fixed by ^Be + p and '^Li + n elastic scattering 
data and a cluster model calculation of Descouvemont and 
Baye using the Volkov II interaction are also shown. The 
calculations have been scaled to best fit the data below 400 
keV. 



well. However, the extrapolated values differ by more 
than 3%. This fact exposes the danger of using high rela- 
tive energy data (above the Ml resonance) to extrapolate 
to zero energy. Here, even limiting the relative energy 
range to < 400 keV and considering only potential models 
that reproduce the available elastic scattering data, the 
extrapolation uncertainty is appreciable. The descrip- 
tion of these data and those considered in the following 
section is significantly better with the potential model 
parameters taken from ^Li -f n rather than ^Be -I- p elas- 
tic scattering data. Since the ^Li -t- n scattering lengths 
are very precisely measured and isospin is expected to be 
a good symmetry in this system, we use the potential- 
model parameters fixed by ^Li -I- n elastic scattering data 
[sij . and regard the difference between the extrapolated 
values using the different potential model parameter sets 
as the extrapolation uncertainty. It is worth noting that 
the ^Be -I- p scattering lengths are so poorly known that 
the ''Li -I- n potential model curve falls well within the la 
upper and lower limit ^Be 4- p scattering length potential 
model curves shown in Fig.|31 Proceeding in this way, we 
obtain 817(0) = 16.6 ± 1.1 (experimental) ± 0.9 (extrap- 
olation) eV b. Several recent measurements have been 
extrapolated to zero energy using a generator coordinate 
method calculation of Descouvemont and Baye that em- 
ployed the Volkov II nucleon-nucleon interaction . As 
Fig. ^ shows, below 400 keV the shapes of this cluster 
model and the ''Be -|- p parameter set potential model are 
extremely similar; the extrapolated zero-energy S factors 
for the 83 MeV/nucleon Coulomb breakup data using 
these two models differ by less than 0.6%. 



IV. CONSISTENCY OF EXPERIMENTAL S17 
DETERMINATIONS 

We fit the data of Filippone et al. , Hammache et 
al. F^, Strieder et al. Jurighans et al. Baby et 
al. and Schiimann et al. [l3| using the same proce- 
dure. The results of these extrapolations of data below 
relative energies of 400 keV using our potential models 
are shown in Table IITTI and the data themselves appear 
in Fig. 13 The Filippone et al. data have been renor- 
malized using the weighted average of their two target 
thickness measurements after adjusting that based on the 
^Li((i, p)®Li reaction using the cross section recommended 
in 19]. They have not been corrected for ®B backscat- 
tering losses. In order to test the mutual consistency of 
the seven radiative capture and Coulomb breakup mea- 
surements considered here, we have treated them all as 
a single data set which we fit with the potential models 
described above. Fitting the 27 data points below 400 
keV in this way, we find values of 46.5 and 46.6 for 26 
degrees of freedom using the ''Li -I- n and central ^Be -f p 
parameter sets respectively. These values correspond 
to a p- value of 0.008, indicating inconsistency within the 
reported uncertainties. Excluding the measurement of 
Junghans et al., which is both very precise and system- 
atically higher than the others, we find values of 11.7 
and 12.6 for 19 degrees of freedom, corresponding to p- 
values of 0.90 and 0.86. This indicates that the other 
six measurements are mutually compatible. The best fit 
potential models for this reduced data set are shown in 
Fig. [SJ and the zero-energy S factors for the ^Li + n and 
^Be -I- p parameter sets are 18.6(4) and 19.2(4) eV b. 
Hence apart from a single discrepant measurement, the 
recent radiative capture and Coulomb breakup measure- 
ments below 400 keV can be described consistently by 
either of our potential model parameter sets. 

Using the potential model parameters taken from the 
1(7 lower and upper limits on the ^Be -f p scattering 
lengths, we obtain 17.6(4) and 19.7(4) eV b for Si7(0) 
from this reduced data set. Regarding the difference be- 
tween the zero energy S factors extrapolated with the dif- 
ferent potential model parameter sets as the la extrapo- 
lation uncertainty and taking the ''Li + n parameter set 
for the central value, we find 817(0) = 18.6 ± 0.4 (ex- 
perimental) ±1.1 (extrapolation) from the six mutually 
consistent radiative capture and Coulomb breakup mea- 
surements. 

Recommended values of 817(0) based on some of the 
same experimental data have been offered in previous 
work 1^ [15I llTj . Our approach differs from others in 
three principal respects. First, our extrapolation proce- 
dure is based on simple potential models constrained by 
elastic scattering data rather than on a complex cluster 
model unconstrained by such data. Second, we limit the 
relative energy range of the experimental data used in 
the extrapolation to < 400 keV. This is important both 
because the shapes of the models differ increasingly as rel- 
ative energy increases, and because the data themselves 
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TABLE IIL Zero energy astrophysical S factors for the ^Be(p, 7)*B reaction obtained by extrapolating radiative capture and 
Coulomb dissociation data at E^ei < 400 keV with potential models constrained by ^Li + n and ^Be + p elastic scattering 
data. The experimental uncertainties are la values. 
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Filippone 


m 


19.1(8) 


19.8(8) 


Strieder 


m 


18.0(14) 


18.7(14) 


Hammache 


[13] 


18.8(8) 


19.5(8) 


Junghans 


[15] 


21.5(3) 


22.2(3) 


Baby 


\M 


19.5(10) 


20.3(10) 


Schumann 




18.4(10) 


19.0(10) 


Davids 


this work 


16.6(11) 


17.1(11) 


All except Junghans 




18.6(4) 


19.2(4) 



are less consistent at high relative energies. Extrapola- 
tions that include higher relative energy data are domi- 
nated by the statistically precise data obtained above the 
Ml resonance, where the uncertainty in the shape of the 
S factor is problematic. Finally, we evaluate the extrapo- 
lation uncertainties by taking the difference between the 
results of potential models with parameters as extreme as 
the precision of the elastic scattering data allow, rather 
than taking the root mean square deviation of a larger 
set of unconstrained models. Although other groups have 
limited their extrapolations to the low relative energy 
regime, we believe that our method has the most solid 
foundation in experimental data, and will yield both the 
most reliable central value and a realistic estimate of the 
extrapolation uncertainty. It is our contention that this 
extrapolation uncertainty has been universally underes- 
timated, and that a more precise measurement of the 
^Be -I- p scattering lengths will be required in order to 
make progress in this area. 

Several values of 817(0) have been deduced from exper- 
imental determinations of the asymptotic normalization 
coefficient (ANC) of the valence proton wave function in 
**B, 17^18) eV b UJ, 17.4(15) eV b and 17.6(17) 
eV b |42|. A weighted average of these three ANC re- 
sults gives 17.4(10) eV b. This is somewhat smaller than 
but perfectly consistent with our central value from the 
radiative capture and Coulomb breakup data. It is in 
particularly excellent agreement with the value deduced 
from the lower limit ''Be + p parameter set potential 
model. 
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FIG. 5: (Color online) Astrophysical S factors for the 
^Be(p, 7)*B reaction from experimental measurements at rel- 
ative energies < 400 keV. The potential models constrained 
by ^Be -|- p and ^Li -I- n elastic scattering data that best fit 
all these data except those of Junghans et al. are also shown. 
See the text for details. 



V. SUMMARY 

In summary, we have studied the reaction model and 
beam energy dependence of the E2 contribution to the 
Coulomb dissociation of '^B through the comparison of 
^Be longitudinal momentum distribution data with first- 
order perturbation theory, continuum-discretized cou- 
pled channels, and dynamical calculations. A coherent 
picture emerges when one considers nine precise direct 
and indirect measurements of the astrophysical S factor 



for the ^Be(p, 7)^B reaction, though one recent measure- 
ment is discrepant with the others. By fitting radiative 
capture and Coulomb breakup data taken below relative 
energies of 400 keV with potential models constrained by 
^Li -|- n and ^Be -I- p elastic scattering data, we examine 
the mutual consistency of recent S17 measurements and 
obtain a recommended value for 817(6) of 18.6 ± 0.4 (ex- 
perimental) ±1.1 (extrapolation) eV b (Icr). This result 
is in good agreement with recent experimental determi- 
nations of the asymptotic normalization coefficient of the 
valence proton wave function in *B. 
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